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Abstract 

In-space docking and separation maneuvers of 
spacecraft that have large fluid mass fractions 
may cause undesirable spacecraft motion in 
response to the impulsive-acceleration-induced 
fluid motion. An example of this potential low- 
gravity fluid management problem arose during the 
£ development of the shuttle/Centaur vehicle. 

Experimentally verified numerical modeling tech- 
niques were developed to establish the propellant 
dynamics, and subsequent vehicle motion, associ- 
ated with the separation of the Centaur vehicle 
from the shuttle orbiter cargo bay. Although the 
shuttle/Centaur development activity has been 
suspended, the numerical modeling techniques are 
available to predict on-orbit liquid motion 
resulting from impulsive accelerations for other 
missions and spacecraft. 

Introduction 

The NASA Lewis Research Center has a long 
history of providing technology for the design 
and operation of fluid systems in the reduced- 
gravity environment of space. * This technology 
program has provided support for the Centaur vehi- 
cle development effort in the form of both analyt- 
ical and experimental studies. Recently the 
Centaur was being modified to make it compatible 
with the shuttle for launch of the Galileo and 
Ulysses spacecraft. During this Centaur redesign 
effort concern arose over the dynamics associated 
with the separation of the Centaur vehicle from 
the shuttle cargo bay. 

Before separation from the shuttle orbiter 
the Centaur vehicle was to be rotated out of the 
cargo bay by a deployment mechanism attached to 
the Centaur integrated support structure (CISS). 
The two vehicles were then to be separated by 
release of several CISS-mounted springs positioned 
circumferential ly around the aft Centaur structure. 
These springs would impart a separation force to 
the Centaur vehicle along its central axis that 
would yield a differential velocity between the 
orbiter and the Centaur of approximately 1 ft/sec. 
However, because of predeployment shuttle maneu- 
vers and atmospheric drag the propellant in the 
Centaur vehicle was not likely to be positioned 
symmetrical ly with respect to the vehicle axis. 
Consequently the resulting spring force at sepa- 
ration would probably not have acted through the 
center of gravity of the Centaur vehicle, and 


rotational motion of the Centaur relative to the 
orbiter could have occurred. A safety concern 
thus arose over whether this possible Centaur 
rotation would be severe enough to cause impact 
between the Centaur and the CISS or the shuttle 
orbiter. 

This paper presents the results obtained from 
an experimental program conducted in the Lewis 
Zero-Gravity Facility to study fluid motion result- 
ing from an impulsive acceleration. The data were 
used to verify numerical modeling techniques that 
were then employed to assess the importance of the 
Centaur separation dynamics concern. The recently 
released NASA-V0F3D computer code, developed at 
the Los Alamos National Laboratory, was employed 
to perform the numerical modeling. This computer 
code provides unique treatment of low-gravity, 
surface-tension-dominated fluid dynamic phenomena 
so that the propellant motion before and after 
the application of the accelerating force can be 
numerically modeled. The computer code predicts 
not only the propellant motion but also the result- 
ing pressure field in the Centaur propellant tanks. 
These computer-generated tank pressure fields can 
then be employed in a vehicle dynamics analysis to 
establish the Centaur motion during separation 
from the shuttle cargo bay. Because of the rela- 
tive weights of the liquid hydrogen and oxygen 
propellants the influence of the hydrogen propel- 
lant motion on the resulting Centaur vehicle 
motion would have been minimal. Consequently the 
experimental and analytical efforts focused on 
predicting the liquid oxygen propellant motion. 

Experimental Program 

The major goal of the experimental program 
was to verify the capability of a numerical anal- 
ysis to predict the behavior of a liquid subjected 
to an impulsive acceleration. An earlier in-house 
experimental program 2 established possible zero- 
gravity locations of the oxygen tank vapor bubble 
before the Centaur separation. The work reported 
on herein examined the subsequent liquid-vapor 
bubble motion resulting from an impulse similar, 
but not scaled, to the separation maneuver. Scal- 
ing was not attempted because the shuttle/Centaur 
program was faced with the challenge of launching 
the Galileo and Ulysses spacecraft during the 
second quarter of 1986. Consequently the program 
used existing hardware so that the experimental 
data could be obtained in time to support the 
projected launch date. 


This paper is declared a work of the U.S. Government and is 
not subject to copyright protection in the United States. 


1 


Copyright 1987 American Institute of Aeronautics and 
Astronautics, Inc. No copyright is asserted in the United States 
under Title 17, U.S. Code. The U.S. Government has a 
royalty-free license to exercise all rights under the copyright 
claimed herein for Governmental purposes. All other rights are 
reserved by the copyright owner. 



Zero-Gravity Facility 

The Zero-Gravity Facility (Fig. 1) consists 
of a concrete-lined 8.7-m-diameter shaft that 
extends 155 m below ground level. A steel vacuum 
chamber, 6.1 m in diameter and 142 m high, is con- 
tained within the concrete shaft. By using the 
Lewis wind tunnel exhaust system in series with 
vacuum pumps located in the facility, the pressure 
in the vacuum chamber is reduced to 13.3 N/m^ 

(1.3x1 atm). The residual air in the drop 
chamber after pumpdown does produce a very low 
drag on the experimental drop package, but this 
force results in an equivalent acceleration acting 
on the package that is estimated to be no greater 
than 10" b g, thus providing an essentially zero- 
gravity environment. Dropping the package from the 
top of the chamber results in approximately 5 sec 
of free-fall time. The package is suspended from 
the top of the vacuum chamber by a support shaft 
(connected to the cylindrical section of the pack- 
age) on a hinged-plate assembly. Once chamber 
pumpdown is completed, the package is released by 
pneumatically shearing a bolt that holds the 
hinged plate in a closed position. Following the 
free-fall period the package is decelerated at the 
bottom of the chamber in a cart filled with small 
pellets of expanded polystyrene. After the drop 
the vacuum chamber is vented to the atmosphere and 
the drop package is returned to ground level. 

Experimental Drop Package 

The experimental drop package (Fig. 2) con- 
sisted of a cylindrical section housing the exper- 
imental apparatus and the electrical systems. A 
closer view of the experimental apparatus is shown 
in Fig. 3 and a schematic diagram in Fig. 4. 

Located below the test tray were 28-V battery 
packs that supplied and regulated direct current 
to the experiment. 

The hardware mounted on the test tray (Fig. 2) 
included a high-speed motion picture camera, a dig- 
ital clock, a lighting system, and the test tank. 
The test tank was mounted on a cradle fastened to 
a linear bearing slide assembly (Fig. 3). During 
the drop a lateral acceleration was applied to the 
test tank by an air cylinder system (Fig. 4). The 
air cylinder consisted of a piston head (inside 
the cylinder walls) and a piston shaft, which 
extended beyond the end of the cylinder. The 
piston shaft was designed to push against the test 
tank frame, which was held in place with a retain- 
ing pin. During the drop the retaining pin was 
retracted by a solenoid. This allowed the prepres- 
surized air cylinder to impart an acceleration to 
the test tank. The time during which the force 
was applied to the tank was controlled by adjust- 
ing the length of the piston stroke. Once the 
piston reached the end of its stroke and no longer 
was in contact with the tank, the tank slid on 
the linear bearing at a constant velocity. 

The length of the linear bearing, and there- 
fore the distance that the tank could travel after 
being accelerated, was dictated by the experimen- 
tal drop package dimensions and the need to view 
the entire carriage displacement. Although it 
would have been desirable, space was not available 
on the slide assembly to mount the camera. How- 
ever, through the use of a wide-angle lens it was 
possible to view approximately 27 cm of container 


motion with the rigidly attached camera. The 
fact that the camera was not mounted on the slide 
assembly made it nearly impossible to precisely 
measure the fluid motion because of both re frac- 
tion and parallax effects (discussed in detail in 
the following section, "Data Analysis"). The tank 
displacement was measured by tracking the motion 
of pointers attached to the test tank, against a 
scale accurate to 1 mm (Fig. 3). The scale was 
mounted above the container in the same plane as 
the tank centerline and in the camera's field of 
view. 

Test Tank and Liquids 

The test tank was an oblate spheroid (ellip- 
soid) formed from a clear acrylic plastic, with 
semimajor and semiminor axes of 2.00 and 1.47 cm, 
respectively (Figs. 5 and 6). This tank geometry 
was chosen for the test program because of its 
similarity with the liquid-oxygen tank used on 
the Centaur vehicle. The tank was fabricated in 
two halves and bolted together with an 0-ring 
seal as illustrated in Fig. 6. 

The liquids used in the zero-gravity tests 
were ethanol and FC-43. Both liquids exhibited a 
near 0° contract angle on the spheroidal tank 
walls. The near 0° contact angle characteristic 
is identical with that of liquid oxygen on the 
Centaur liquid-oxygen-tank surface. Other key 
properties of the test liquids and liquid oxygen 
(i.e., surface tension, density, and viscosity) 
are presented in Table I. 

It would have desirable to select fluids for 
the experimental program that had viscosity close 
to that of liquid oxygen. However, the results 
presented in Ref. 2 indicate that low-viscosity 
fluids take so long to form the zero-gravity 
liquid-vapor interface configuration that the 
5-sec drop time would not be sufficient to allow 
the subsequent acceleration of the test tank. 
Consequently test fluids with a significant vari- 
ation in viscosity were selected to determine 
whether the numerical analysis correctly accounted 
for viscosity effects. 

Test Setup and Procedure 

Before each series of tests the test tank was 
cleaned ultrasonical ly so that surface contamin- 
ation was avoided. Immediately before each test 
the tank was rinsed with a solution of distilled 
water, dried in a warm air dryer, and then rinsed 
with the test liquid. After rinsing, the tank was 
filled to the desired liquid level with a syringe, 
hermetically sealed, and mounted on the test 
vehicle. 

The entire test vehicle was balanced and the 
cradle on which the test tank was mounted was 
leveled to ensure that the ellipsoidal tank was 
oriented properly before and during the drop. 

After the tank was mounted and the vehicle bal- 
anced, the entire drop package was suspended from 
the top of the facility. The facility vacuum 
chamber was then closed and pumped down and the 
drop package was released. At package release the 
lights and the motion picture camera were activated. 
Approximately 2 sec after release lateral acceler- 
ation was applied to the experimental apparatus, 
and the test tank traversed the linear bearing. 
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Data Analysis 

The data obtained in this investigation were 
collected daring the free-fall period by a high- 
speed motion picture camera. Information on the 
zero-gravity liquid-vapor interface configuration 
was taken from the film by using a film analyzer. 
From the analyzer the observed interface shape 
could be directly plotted at various times through- 
out the drop. In the test films the general out- 
line of the interface was defined by a dark band. 
The points used for plotting the interface shape 
followed the outer perimeter of this dark band 
(Fig. 6). Because of symmetry of the interface 
about the ellipsoid centerline before and during 
the drop, this outer perimeter corresponded to 
the part of the interface located in the plane 
defined by the ellipsoid centerline and the tank 
minor axis. 

This observed interface, however, became dis- 
torted by refraction effects and at times because 
of the angle at which the camera viewed the test 
tank. The constantly present refraction of the 
observed interface was caused by the walls of the 
ellipsoidal tank and a layer of test liquid. 
Another form of distortion was created when the 
container was moved during a drop. When this 
occurred, the tank started at one end of the 
linear bearing track and traveled to the other 
end. The camera did not move with the tank 
(because of hardware limitations discussed pre- 
viously) but instead was mounted approximately 
halfway along the distance the test tank traveled. 
This caused a parallax effect when the tank was 
not at the center of the track. It was possible 
to correct for the refraction error alone by using 
a calibration grid from which actual interface 
coordinates could be determined (see Ref. 3 for a 
detailed explanation). Correcting for parallax 
effects in combination with refraction, however, 
was not possible. Therefore, when the test tank 
was near the ends of the track, an approximate, or 
qualitative, liquid-vapor interface (location and 
configuration) was determined by studying the test 
film and using the calibration grid. When the 
tank was near the middle of the track and was 
positioned along the optical axis of the camera 
lens, the calibration grid could be employed to 
obtain an actual interface configuration. 

In addition to the liquid-vapor interface 
plot the lateral acceleration characteristics of 
the tank were also obtained from the test film. 
This was accomplished by recording lateral dis- 
placement versus time data from the film. Two 
pointers (Fig. 3) were attached to the tank and 
used for reading displacement measurements. One 
pointer was aligned with the camera centerline 
when the tank was at the far end of the track 
before the application of the impulsive acceler- 
ation. The other pointer became centered in the 
camera's field of view as the tank traveled across 
the linear bearing track. The use of two pointers 
and a correction for parallax effects allowed for 
accurate displacement readings regardless of the 
location of the tank. 

Accurate knowledge of the time the tank first 
started to accelerate was critical in correlating 
the initial displacement characteristics with time. 
In this investigation the starting time was indi- 
cated by illumination of a light-emitting diode. 
This light-emitting diode was activated when the 


retaining pin was rotated away from the tank base, 
allowing the air piston to accelerate the tank. 
Since the stroke of the piston was usually 
adjusted to only 1.27 cm, the tank experienced a 
lateral acceleration for 0.1 to 0.2 sec. 

The magnitude of the acceleration was deter- 
mined by plotting displacement (OISP) versus time 
points on log-log graph paper. A curve was fitted 
to these points by using a slope of 2; the accel- 
eration (ACC) was then obtained from the appro- 
priate intercept of the curve (Fig. 7 is a typical 
plot). After the piston reached the end of its 
stroke, the acceleration value went to zero and 
the velocity of the test tank remained nearly con- 
stant. This velocity was determined by plotting 
displacement versus time and taking the slope of 
the plot (Fig. 8) . 

The variation of experimental conditions 
encompassed by the five tests conducted as part of 
this program and the anticipated shuttle/Centaur 
deployment conditions are shown in Table II. Also 
included in Table II are the calculated values of 
the test tank accelerations and the total time of 
the accelerating impulses. 

Test Results and Discussion 

In order to understand the problems associ- 
ated with the dynamics of low-gravity liquid- 
vapor systems in ellipsoidal tanks, the initial 
liquid-vapor interface location before any dis- 
turbance must be known. The initial phase of the 
test program^ was performed with the objective of 
determining equilibrium zero-gravity interface con- 
figurations that could exist in an ellipsoid. 

During this first phase of the study the test tank 
was not given an acceleration but instead was 
clamped to the linear bearing in a fixed position 
such that the tank was aligned with the optical 
axis of the motion picture camera (Fig. 2). The 
zero-gravity environment (drag on the experimental 
package has an insignificant effect on the 
results), coupled with solid-liquid-vapor combina- 
tions that produced near 0° static contact angles, 
resulted in spherical liquid-vapor interface shapes 
inside the ellipsoidal tank. The experimental 
work that established equilibrium zero-gravity 
liquid-vapor interface shapes and locations is 
discussed in detail in Ref. 2. 

The results described in Ref. 2 indicate 
that, depending on the size of the vapor bubble, 
the liquid-vapor interface can exist at several 
locations in an ellipsoidal tank. This is true as 
long as the interface is able to meet the minimum 
energy requirement by forming a sphere within the 
walls of the tank. For the shuttle/Centaur sepa- 
ration dynamics problem, only liquid fill levels 
greater than 70 percent were of interest. Liquid 
fill levels of 80 percent or more resulted in vapor 
bubbles with diameters small enough so that they 
could be accommodated anywhere within the tank 
boundaries. As the interface radius of curvature 
increased with reduced fill levels, between 80 and 
70 percent, the vapor bubble tended to form in the 
center of the ellipsoidal tank since this was the 
only location where a spherical bubble could fit 
within the tank boundaries. 

After a better understanding of the equilib- 
rium zero-gravity interface configurations had 
been acquired, 2 the experimental investigation 
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continued by applying a lateral acceleration to 
the ellipsoidal tank. The orientation of the tank 
relative to the gravitational vector (Fig. 9) was 
selected so that the resulting zero-gravity vapor 
bubble would be at the assumed worst-case location 
(i.e., at the farthest extreme from the center of 
the tank, thus providing a maximum center-of- 
gravity offset) before the start of the impulsive 
acceleration for tests 1, 2, and 3. Test condi- 
tions, liquids, and fill levels were chosen so 
that the time to form a nearly quiescent interface 
was as short as possible; usually less than 2 sec. 
During the test sequence the start of the lateral 
acceleration was delayed by 2.0 to 2.5 sec to 
allow the liquid-vapor interface to form its zero 
Bond number configuration. This formation period 
was sufficiently short to allow adequate time to 
view the lateral acceleration of the test tank and 
the resulting motion of the tank and the test 
fluid. 

Figures 9 and 10 show a sequence of pictures 
from one of the zero-gravity tests (test 1, 

Table II) in which the ellipsoidal tank was accel- 
erated. In Fig. 9 the first portion of the drop 
can be seen. This part of the drop involved the 
formation of a stable liquid-vapor interface. The 
initial condition of the experimental drop package 
was such that the liquid bottomed in the tank and 
was symmetric about the container major axis 
(Fig. 9(a)). Upon entry into the low-gravity 
environment the liquid-vapor interface underwent 
an oscillatory motion but remained symmetric about 
the ellipsoid major axis (Fig. 9(b)). This tran- 
sient motion displaced the bubble along the major 
axis (Fig. 9(c)) to the other side of the tank, 
where the liquid motion damped and the vapor bub- 
ble became stable (Fig. 9(d)). The process just 
described was well understood from the first phase 
of the experimental work and took approximately 
2.5 sec. 

Once the interface became quiescent, the 
lateral acceleration was applied to the test tank. 
Just after the acceleration was applied, the vapor 
bubble started to move from the center of the tank 
toward the left wall (Fig. 10(a)). (The force was 
applied to the right side of the tank as viewed 
from the camera.) After the force had been applied 
for approximately 0.1 sec, a geyser appeared to 
form and the vapor bubble became highly deformed 
(Fig. 10(b). When the piston reached the end of 
its stroke and the applied force was no longer 
present, the geyser collapsed - leaving an inter- 
face that appeared to have no distinct shape 
(Fig. 10(c)). However, in a very short time 
(approximately 0.4 sec) after the acceleration 
became zero, the interface again appeared to be 
dominated by capillary forces and started to take 
on a spherical shape (Fig. 10(d)). As the ellip- 
soidal tank continued to move across the track, 
the interface traveled from the left side to the 
right side of the tank (Fig. 10(e)). Toward the 
end of the drop period the liquid motion damped, 
and the liquid-vapor interface appeared to be 
forming its equilibrium spherical zero-gravity 
configuration (Fig. 10(f)). 

Figure 11 shows plots of the liquid-vapor 
interface corresponding to some of the photographs 
just described. Figure 11(a) illustrates the 


initial condition of the low-gravity interface 
before any impulse. This plot was derived from 
the first phase of the experimental program. 

Figures 11(b) and (d) correspond to times when 
the interface was distorted due to refraction and 
parallax and therefore are qualitative interpre- 
tations of the observed interface (determined as 
described in the section "Data Analysis*'). A more 
accurate interface shape was plotted when the tank 
was near the middle of the track (Fig. 11(c)). 

Note that this plot is not complete, however, 
because refraction was particularly high near the 
edge of the ellipsoidal tank. 

Four other impulsive acceleration runs were 
performed and each run had similar results. 

Figures 12 to 15 are plots, similar to Fig. 11, of 
the liquid-vapor interface configuration at key 
times during test runs 2 to 5, respectively. In 
every run but one (test 3, which had the lowest 
impulsive acceleration) a geyser formed, and in 
all tests the vapor bubble remained intact and was 
becoming stable near the end of the test period. 

This part of the program was performed to 
provide data for the verification of a computer 
code and did not actually scale Centaur tank con- 
figuration or flight conditions. Time constraints 
permitted only a limited number of tests to be 
made. Consequently no attempt was made to develop 
correlating parameters relating the influence of 
the impulsive lateral acceleration to the observed 
behavior of the liquid and the vapor within the 
test tank. 

Numerical Modeling 

The NASA-V0F3D computer code^ was used to 
provide (1) an analytical comparison with the 
results from three of the previously described 
tests, thus providing a variation in both liquid 
fill level and fluid properties, and (2) an analyt- 
ical prediction of the Centaur liquid oxygen motion 
resulting from the deployment of the vehicle from 
the shuttle cargo bay, assuming a worst-case ini- 
tial location for the vapor bubble. The NASA-V0F3D 
program simulates three-dimensional incompressible 
flows with free surfaces by using the volume-of- 
fluid algorithm. It is specifically designed to 
predict liquid flows in a low-gravity environment, 
where the free surface physics must be accurately 
treated. The computational technique is based on 
the use of donor-acceptor differencing to track 
the free surface across an Eulerian grid. 

NASA-V0F3D is the latest of a family of com- 
puter codes developed at the Los Alamos National 
Laboratory under an interagency agreement with 
NASA Lewis*. The first of these codes, NASA 
S0LA-V0F,^ is a two-dimensional program written 
specifically for simulating the axisymmetric drain- 
ing of spacecraft propellant tanks. The recently 
released NASA-V0F2D 6 and NASA-V0F3D programs con- 
tain several improvements over the NASA S0LA-V0F 
code. 


*A11 three computer codes are available, for 
a nominal fee, from the National Energy Software 
Center, Argonne National Laboratory, Argonne, 
Illinois. 
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The NASA-V0F3D program, like its NASA-VOF2D 
predecessor, is highly structured so that indi- 
vidual components may be easily modified to fit 
specific problems or to accept subsequent code 
upgrades. This approach allows greater efficiency 
in operation, simplifies problem setup, and facil- 
itates code modification. Although NASA-V0F3D can 
run many problems without alteration, a specific 
application will generally require program alter- 
ations for special inflow or outflow ports, com- 
plicated geometries, or unusual initial conditions. 
Seldom does anything have to be deleted; rather 
the modifications are usually additions to the 
existing structure. 

Because much of the numerical methodology is 
very similar in two and three dimensions, the reader 
will frequently be referred to the NASA-V0F2D 
report 6 for a fuller description of topics men- 
tioned in this paper. In particular. Ref. 6 dis- 
cusses the developments leading to the use of 
SOLA-VOF methodology in low-gravity liquid flows. 
NASA-V0F3D incorporates the NASA-V0F2D improve- 
ments in donor-acceptor differencing, in the abil- 
ity to accurately model curved boundaries, and in 
the inclusion of a conjugate residual option as 
well as a successive-overrelaxation (SOR) option 
for advancing the pressures (and velocities) in 
time throughout the computing mesh. The conjugate 
residual method is vectorized for the Cray com- 
puter and employs a scaled coefficient matrix. It 
was found to be necessary for accuracy to radi- 
cally alter the procedure for calculating three- 
dimensional surface tension effects, as will be 
explained in the following section. 

The basic NASA-VOF algorithm is based on the 
solution of the Navier-Stokes equations for an 
incompressible fluid. These equations apply at 
every point inside the fluid. The complete 
Navier-Stokes equations for an incompressible 
fluid are solved by finite differences, with the 
surface tension effects included. Several sets 
of quantities must be incremented to advance the 
calculation through one time step. One set of 
procedures is applied to all computational cells 
open to flow; other procedures are applied to 
fluid cells only, with alternative procedures for 
surface cells. 

The free surfaces are treated by introducing 
a function F that is defined to be unity at any 
point occupied by the liquid, and zero elsewhere. 
When averaged over a cell of the computing mesh, 
the value of F represents the fractional volume 
of the cell occupied by liquid. Cells with F 
values between one and zero are partially filled 
with liquid; they are either intersected by a free 
surface or contain voids (bubbles) smaller than 
cell mesh dimensions. The F function is used to 
determine which cells contain a free surface bound- 
ary and where the liquid is located in the cells. 
Gradients of F determine the mean local surface 
normal and together with F values permit con- 
struction of an approximate interface. Surface 
tension and wall adhesion forces may then be com- 
puted. The special requirements for construction 
of gradients of the step function F are dis- 
cussed in Ref. 6. 

Numerical Methodology 

This section briefly describes the numerical 
methods employed by the NASA-VOF computer codes. 


Only the key features and limitations of the pro- 
grams are presented; for more detail examine 
Refs. 4 and 6. The programs use an Eulerian frame 
and allow for an arbitrary number of free surface 
segments with any reasonable shape. Variable mesh 
spacing is permitted in all coordinate directions, 
as is necessary for efficient numerical solution 
of many problems. However, the variable mesh capa- 
bility is subject to reasonable requirements for 
avoiding neighboring cells of grossly different 
size and cells of large aspect ratio. 

Typical boundary conditions for the computa- 
tional mesh are imposed in a standard way by the 
use of flags denoting the tank walls. Other types 
of boundary conditions, if needed, must be sup- 
plied by the user and placed in a special subrou- 
tine. An example of the use of special boundary 
conditions (i.e., a tank drain hole) is contained 
in the NASA SOLA-VOF report. 

The program specifies the location of obsta- 
cles and the initial locations of liquid, vapor, 
and liquid-vapor interfaces. When appropriate, a 
liquid-vapor interface may be calculated as an 
equilibrium surface; otherwise its shape must be 
specified by the user in advance. 

The initial liquid velocities are specified 
as constant or zero. It would, in principle, be 
possible to take the liquid to be initially in 
nonuniform motion. However, it is very rare for 
a dynamic velocity field to be known with suffi- 
cient precision to be suitable as an initial con- 
dition. Although the pressure field is initially 
zeroed, this command may be overridden to produce 
a hydrostatic head as the system's initial pres- 
sure field. It is important that an accurate 
initial pressure field be employed since problems 
with poor initial pressure definition may have 
difficulty. 

All the surface physics algorithms developed 
for NASA-V0F2D were improved and generalized dur- 
ing the development of the three-dimensional 
modeling capability. In two dimensions one sur- 
face curvature can be calculated algebraically, 
and one curvature must be computed by finite- 
difference approximations. Moreover the algebra- 
ically determined curvature is typically the 
larger one in the crucial near-axis cells. In 
three dimensions both surface curvatures must be 
computed by less accurate finite-difference 
approximations. In addition, considerable refine- 
ment of the two-dimensional surface tension algo- 
rithms was necessary during the development of 
the three-dimensional computer code. Reference 4 
contains a detailed discussion of how these sur- 
face physics effects were accounted for 
numerically. 

Considerable effort was expended in obtaining 
a realistic wall adhesion algorithm for two- 
dimensional problems (it is a standard feature of 
NASA-V0F2D) . This turned out to be of crucial 
importance in calculating such flow features as 
the leading-edge liquid velocity in a propellant 
reorientation problem. In two dimensions the task 
was accomplished (1) by defining a cell index 
array that specifies whether the wall adhesion 
option or the surface tension option is to be 
selected and (2) by defining angles, consistent 
with the tank geometry, that determine the direc- 
tion in which any wall adhesion force is to act. 


5 


In two dimensions the wall adhesion calculation is 
simplified because the force points either up or 
down the tank wall. This approach neglects the 
effect of liquid contact angle since the contact 
angle is usually very small. In contrast, a suc- 
cessful three-dimensional algorithm must introduce 
an additional angle, redefined for each computa- 
tional cycle, that determines in what direction 
the wall-adhesion force is to point. 

The present version of the three-dimensional 
program does not permit the accurate calculation 
of wall adhesion except for straight-walled cells; 
in general, it forces the cell indicator flags to 
choose the surface tension option. The inclusion 
of a realistic wall adhesion calculation into the 
surface physics routines is one of the principal 
tasks that must be accomplished if NASA-V0F3D is 
to evolve into a truly general-purpose code. How- 
ever, it is emphasized that the present program 
still has a very significant range of applicability 
that permits a wide range of design questions to 
be meaningfully addressed. For example, questions 
concerning the bulk motion of fluid in a given 
geometry may be addressed without the interaction 
of the liquid and its tank wall being predicted 
with high accuracy. 

Analytical Results and Discussion 

Three of the Zero-Gravity Facility tests 
(tests 1, 4, and 5 in Table II) were numerically 
modeled. Inputs to the NASA-V0F3D computer code 
included liquid properties (Table I), tank geom- 
etry, vapor bubble size and location before the 
start of impulsive acceleration, and the duration, 
direction, and magnitude of the applied 
acceleration. 

The results of the numerical modeling activ- 
ity for test 1 are shown in Fig. 16 as cross- 
sectional views through the center of the tank. 

The two-dimensional views of the vapor bubble and 
the liquid velocity vectors shown are a graphical 
output feature of the NASA-V0F3D computer code. 

For each time step the liquid velocity vectors are 
scaled to the maximum liquid velocity at that 
time. The sequence in time of the first, third, 
and last cross-sectional views of the vapor bubble 
shown in Fig. 16 were selected to correspond 
approximately in time with the plots of the exper- 
imental data presented in Figs. 11(a) to (c). The 
numerical modeling was terminated at approximately 
0.6 sec (Fig. 16(f)) when significant vapor bubble 
motion was no longer observed and the vapor bubble 
began to break up. 

The numerical modeling results corresponding 
to tests 4 and 5 are presented in Figs. 17 and 18, 
respectively. Once again, the sequence in time of 
the cross-sectional views of the vapor bubbles was 
chosen to correspond approximately in time with 
the plots of the liquid-vapor interfaces derived 
from the experimental data (Figs. 14 and 15). 

Comparing the experimental data generated in 
the Lewis Zero-Gravity Facility with the NASA- 
V0F3D numerical modeling output led to the con- 
clusion that the computer code could be used to 
predict liquid-oxygen motion in the Centaur vehi- 
cle. Excellent agreement existed between the 
experimental and numerical results for the dur- 
ation of the accelerating force. During this time 


most of the liquid motion took place and the high- 
est liquid velocities were calculated. At later 
times, in contrast to the experimental observa- 
tions, the numerically generated ullage fragmented, 
probably reflecting the continuing difficulty in 
accurately modeling surface-tension-dominated fluid 
dynamic phenomena. However, before the ullage 
broke up, the liquid motion had significantly sub- 
sided, as observed experimentally and determined 
numerically. Since the liquid momentum flux was 
proportional to the square of the liquid velocity, 
this numerical inaccuracy should not affect the 
ultimate objective, the prediction of the liquid- 
oxygen propellant and Centaur vehicle motion in 
response to an impulsive acceleration. 

Similar to the numerical modeling of the 
experimental data from the Zero-Gravity Facility, 
the NASA-V0F3D computer code was employed to pre- 
dict the liquid oxygen motion in the Centaur 
vehicle' during a deployment maneuver. Inputs 
to the code included liquid-oxygen properties 
(Table I), the vapor bubble geometry and its 
assumed worst-case initial location (maximum dis- 
placement of the liquid center of gravity from 
the Centaur vehicle centerline), the tank geometry, 
and the characteristics of the applied acceler- 
ation (Table II). 

The results of the NASA-V0F3D modeling effort 
for the shuttle/Centaur program are shown in 
Figs. 19 to 23. The left portion of each figure 
is a side view of the liquid-oxygen tank showing a 
cross section through the Centaur vehicle center- 
line. (Note that, in contrast with the oblate 
spheroid tank configuration used for the experi- 
mental tests, the shuttle/Centaur liquid-oxygen 
tank has been enlarged by the addition of a cylin- 
drical midsection so that it appears to be nearly 
spherical.) The right side of each figure is a 
cross-sectional view of the tank that is perpen- 
dicular to the Centaur vehicle centerline and 
taken in a plane approximately one-third of the 
tank height from the top of the tank. Figure 19 
shows the assumed worst-case initial location for 
the vapor bubble. One second after the start of 
the 0.44-sec impulsive acceleration the liquid 
velocity reached a maximum of approximately 
50 cm/sec (Fig. 20). As time progressed (Figs. 21 
to 23) the liquid velocities decayed and the vapor 
bubble became highly distorted, nearly separating 
into two distinct regions. 

The numerical output of the NASA-V0F3D com- 
puter code included the prediction of the pressure 
distribution in the Centaur liquid-oxygen tank as 
well as the liquid velocity field shown graphically 
in Figs. 20 to 23. These pressure distribution 
predictions can be employed in a vehicle dynamics 
analysis to establish the relative motion of the 
Centaur and shuttle vehicles during a deployment 
maneuver. For the first two planned Centaur-in- 
shuttle missions, the launches of the Galileo and 
Ulysses spacecraft, no hazardous vehicle motion 
was predicted. 

Concluding Remarks 

The NASA Lewis Research Center, during the 
past two decades, has been a leading developer of 
low-gravity fluid management technology. Recent 
analytical and experimental efforts have addressed 
the potential problem associated with the separa- 
tion of the Centaur vehicle from the shuttle 


6 



orbiter. Experimental data obtained in the unique 
Lewis Zero-Gravity Facility were used to assess 
the capabilities of the NASA-V0F3D computer code, 
developed under an interagency agreement with the 
Los Alamos National Laboratory. The NASA-V0F3D 
program was then employed to predict the liquid- 
oxygen propellant motion and pressure field within 
the Centaur tankage during the planned deployment 
maneuver. 

Based on the experimental and analytical work 
performed it was concluded that no hazardous vehi- 
cle motion would have resulted from the Centaur 
separation maneuver if the propellant tanks were 
nearly full. Although the shuttle/Centaur devel- 
opment activity has been suspended, the problem 
solution presented serves to illustrate the capa- 
bilities of the NASA-V0F3D computer code and the 
ability of Lewis to provide timely fluid manage- 
ment technology information to developers of 
advanced space systems. 
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TABLE I. - LIQUID PROPERTIES 


Liquid 

Surface 

tension, 

dyn/cm 

Density, 

g/cm J 

Dynami c 
viscosity, 
cP 

Ethanol 
(at 20 °C) 

22.3 

0.79 

1.20 

FC-43 

(at 20.5 °C) 

16.8 

1.90 

6.41 

Liquid oxygen 
(at -180 *C) 

10.0 

1.07 

.19 


TABLE II. - SUMMARY OF ZERO-GRAVITY FACILITY EXPERIMENTAL 
CONDITIONS AND ANTICIPATED SHUTTLE/CENTAUR DEPLOYMENT 
CONDITIONS FOR GALILEO MISSION 


Test 

Liquid 

Fill 
level , 
percent 

Applied 

acceleration, 

9 

Time 

acceleration 

applied, 

sec 

1 

Ethanol 

90 

0.25 


2 


90 

.22 


3 


90 

.10 


4 


75 

.22 


5 

FC-43 

90 

.14 

.132 

Shuttle/Centaur 

Liquid 

oxygen 

89 

40.133 

0.440 


a Maximum value (applied acceleration = 0.133 cos 3.47 t, where 
t denotes time in seconds.) 
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•Concrete-lined shaft 
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Bottom of shaft 155 meters 
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FIGURE 3. - EXPERIMENTAL APPARATUS MOUNTED ON TEST TRAY. 
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FIGURE 1. - NASA LEWIS RESEARCH CENTER ZERO GRAVITY FACILITY. 


.'•s— 6. 1-Meter -diameter steel 
vacuum chamber 


FIGURE 2. - ZERO-GRAVITY FACILITY EXPERIMENTAL DROP PACKAGE. 
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original page fS 

OE POOR QUALITY 


TIMER AND 



GENERATOR 

FIGURE - SCHEMATIC OF EXPERIMENTAL APPARATUS. 


DIRECTION 

EXPERIMENTAL 

PACKAGE 

DROPPED 

ELLIPSOID CENTERLINE 7 



FRONT VIEW RELATIVE TO SIDE VIEW RELATIVE TO 

CAMERA CAMERA 

FIGURE 5. - ELLIPSOIDAL CONFIGURATION OF MODEL TANK. 



FIGURE 6. - EQUILIBRIUM LI QUID- VAPOR INTERFACE CONFIGURATION IN 
ELLIPSOIDAL TANK. BOND NUMBER, 0; FILL LEVEL, 80 PERCENT. 



TIME, SEC 


FIGURE 7. - REPRESENTATIVE LATERAL ACCELERATION CHAR- 
ACTERISTICS USING ASSUMED SLOPE OF 2. LATERAL AC- 
CELERATION FOR THIS PLOT, 0.22 g. 



FIGURE 8. - REPRESENTATIVE LATERAL VELOC- 
ITY CHARACTERISTIC. LATERAL VELOCITY 
FOR THIS PLOT, 22.9 CM/SEC. 
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OF POOR QUALITY 


(A) INITIAL LIQUID-VAPOR INTERFACE CONFIGURATION; TIME, 0 SEC 


(B) TRANSIENT INTERFACE SHAPE; TIME, 0.09 SEC INTO DROP 


(C) DISPLACEMENT OF INTERFACE ALONG TANK MAJOR AXIS; TIME, 0.39 
SEC INTO DROP. 


(D) EQUILIBRIUM ZERO-GRAVITY INTERFACE CONFIGURATION; TIME 
2.50 sec INTO DROP. 


FIGURE 9. - TRANSIENT MOTION OF INTERFACE UPON ENTERING LOW GRAVITY. FILL LEVEL, 90 PERCENT. | (TANK MAJOR AXIS PARALLEL TO DROP DIRECTION.) 



ORIGINAL PAGE 13 

OF POOR QUALITY 



(E) MOVEMENT OF INTERFACE ALONG TANK MAJOR AXIS DURING CONSTANT- (F) INTERFACE AGAIN DOMINATED BY CAPILLARY FORCES; TIME, 1.09 SEC. 

VELOCITY PERIOD OF TANK MOTION; TIME, 0.62 SEC. 


FIGURE 10. - LIQUID-VAPOR INTERFACE CONFIGURATION DURING AND AFTER IMPULSIVE ACCELERATION. TEST 1; FILL LEVEL, 90 PERCENT; IMPULSE, 0.25g 
FOR 0.101 SEC. (TIME MEASURED FROM START OF IMPULSIVE ACCELERATION.) 
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(C) QUANTITATIVE VIEW OF LIQUID-VAPOR INTERFACE: 
TIME, 0.62 SEC. 



(B) QUALITATIVE VIEW OF LIQUID-VAPOR INTER- 
FACE; TIME, 0.10 SEC. 



(D) QUALITATIVE VIEW OF LIQUID-VAPOR INTER- 
FACE; TIME, 1.09 SEC. 


FIGURE 11. - PLOTS OF LIQUID-VAPOR INTERFACE CONFIGURATION - TEST 1. FILL LEVEL, 90 PERCENT: IMPULSE, 0.25 g FOR 0.101 sec. 
(TIME MEASURED FROM START OF IMPULSIVE ACCELERATION.) 




FORCE 

APPLIED 




(B) QUALITATIVE VIEW OF LIQUID-VAPOR INTERFACE; 



(C) QUANTITATIVE VIEW OF LI QUID- VAPOR INTERFACE: 
TIME, 0.69 SEC. 


(D) QUALITATIVE VIEW OF LIQUID-VAPOR INTER- 
FACE; TIME, 1.18 sec. 


FIGURE 12. - PLOTS OF LIQUID-VAPOR INTERFACE CONFIGURATION - TEST 2. FILL LEVEL, 90 PERCENT; IMPULSE, 0.22 g FOR 0.108 sec. 
(TIME MEASURED FROM START OF IMPULSIVE ACCELERATION.) 
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(B) QUALITATIVE VIEW OF LIQUID-VAPOR INTERFACE; 
TIME. 0.15 sec. 



(C) QUANTITATIVE VIEW OF LIQUID-VAPOR INTERFACE; (D) QUALITATIVE VIEW OF LIQUID-VAPOR INTERFACE; 

TIME. 0.93 sec. TIME. 1.72 sec. 

FIGURE 13. - PLOTS OF LIQUID-VAPOR INTERFACE CONI .URATION - TEST 3. FILL LEVEL. 90 PERCENT; IMPULSE, 0.10 g FOR 0.160 sec. 
(TIME MEASURED FROM START OF IMPULSIVE ACCELERATION.) 
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(A) EQUILIBRIUM ZERO-GRAVITY CONFIGURATION. (B) QUALITATIVE VIEW OF LIQUID-VAPOR INTERFACE; 

TIME, 0.10 SEC. 



(C) QUANTITATIVE VIEW OF LIQUID-VAPOR INTERFACE; (D) QUALITATIVE VIEW OF LIQUID-VAPOR INTERFACE; 

TIME, 0.67 SEC. TIME, 1.18 sec. 


FIGURE 14. - PLOTS OF LIQUID-VAPOR INTERFACE CONFIGURATIONS - TEST 4. FILL LEVEL, 75 PERCENT; IMPULSE, 0.22 g FOR 0.109 SEC. 
(TIME MEASURED FROM START OF IMPULSIVE ACCELERATION.) 









(A) LIQUID-VAPOR INTERFACE CONFIGURATION (B) INTERFACE CONFIGURATION AND LIQUID VE- <C) INTERFACE CONFIGURATION AND LIQUID VE- 

BEFORE START OF IMPULSIVE ACCELERATION. LOCITY VECTORS 0.05 sec AFTER IMPULSE LOCITY VECTORS 0.10 SEC AFTER IMPULSE IN- 
INITIATION (MAXIMUM VELOCITY **15 cm/sec). ITIATION (MAXIMUM VELOCITY = 27 CK/SEC). 



(D) INTERFACE CONFIGURATION AND LIQUID (E) INTERFACE CONFIGURATION AND LIQUID (F) INTERFACE CONFIGURATION AND LIQUID VE- 

VELOCITY VECTORS 0.25 SEC AFTER IMPULSE VELOCITY VECTORS 0.35 SEC AFTER IMPULSE LOCITY VECTORS 0.60 SEC AFTER IMPULSE INI- 

INITIATION (MAXIMUM VELOCITY = 10 CM/SEC). INITIATION (MAXIMUM VELOCITY * 7 CM/SEC). T I AT I ON (MAXIMUM VELOCITY * 5 CM/SEC). 


FIGURE 16. - NUMERICAL MODELING OF LIQUID AND VAPOR MOTION IN RESPONSE TO IMPULSIVE ACCELERATION - TEST 1. FILL LEVEL, 90 PERCENT; IMPULSE, 
0.25 g FOR 0. 101 sec. 
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(A) LIQUID-VAPOR INTERFACE CONFIGURATION (B) INTERFACE CONFIGURATION AND LIQUID VE- (C) INTERFACE CONFIGURATION AND LIQUID 

BEFORE START OF IMPULSIVE ACCELERATION. LOCITY VECTORS 0.05 SEC AFTER IMPULSE VELOCITY VECTORS 0.10 SEC AFTER IMPULSE 

INITIATION (MAXIMUM VEL0CITY=16 CM/SEC). INITIATION (MAXIMUM VELOCITY =27 cm/sec). 



(D) INTERFACE CONFIGURATION AND LIQUID VE- (E) INTERFACE CONFIGURATION AND LIQUID VE- (F) INTERFACE CONFIGURATION 0.65 SEC AFTER 

LOCITY VECTORS 0.15 SEC AFTER IMPULSE LOCITY VECTORS 0.35 SEC AFTER IMPULSE IMPULSE INITIATION. BULK LIQUID VELOCITY 

INITIATION (MAXIMUM VELOCITY =14 CM/SEC). INITIATION (MAXIMUM VEL0CITY = 5 CM/SEC). VECTORS NEARLY ZERO. 


FIGURE 17. - NUMERICAL MODELING OF LIQUID AND VAPOR MOTION IN RESPONSE TO IMPULSIVE ACCELERATION - TEST 4. FILL LEVEL. 75 PERCENT; IMPULSE, 
0.22 g FOR 0.104 SEC. 
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(A) LIQUID-VAPOR INTERFACE CONFIGURATION 
BEFORE START OF IMPULSIVE ACCELERATION. 



(B) INTERFACE CONFIGURATION AND LIQUID VE- 
LOCITY VECTORS 0.05 SEC AFTER IMPULSE 
INITIATION (MAXIMUM VELOCITY* 10 cm/sec). 


(C) INTERFACE CONFIGURATION AND LIQUID VE- 
LOCITY VECTORS 0.13 SEC AFTER IMPULSE 
INITIATION (MAXIMUM VELOCITY* 22 CM/SEC). 




(D) INTERFACE CONFIGURATION AND LIQUID VE- 
LOCITY VECTORS 0.25 sec AFTER IMPULSE IN- 
ITIATION (MAXIMUM VELOCITY-7 cm/sec). 



(E) INTERFACE CONFIGURATION AND LIQUID VE- 
LOCITY VECTORS 0.45 SEC AFTER IMPULSE 
INITIATION (MAXIMUM VEL0CITY = 6 CH/SEC). 



(F) INTERFACE CONFIGURATION AND LIQUID VE- 
LOCITY VECTORS 0.85 SEC AFTER IMPULSE 
INITIATION (MAXIMUM VELOCITY =4 CM/SEC). 


FIGURE 18. - NUMERICAL MODELING OF LIQUID AND VAPOR MOTION IN RESPONSE TO IMPULSIVE ACCELERATION - TEST 5. FILL LEVEL, 90 PERCENT; IMPULSE, 
0.14 g FOR 0.132 sec. 
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(A) CROSS-SECTIONAL VIEW OF LIQUID-OXYGEN TANK. (B) SECTION A-A. 


FIGURE 19. - NUMERICAL MODELING OF CENTAUR LIQUID-OXYGEN PROPELLANT MOTION DURING SEPARATION MANEUVER. INITIAL 
ASSUMED WORST-CASE LIQUID-VAPOR INTERFACE CONFIGURATION; FILL LEVEL, 89 PERCENT, IMPULSE APPLIED FOR 0.44 sec 
WITH MAXIMUM ACCELERATION OF 0.133 g. 



(A) CROSS-SECTIONAL VIEW OF LIQUID-OXYGEN TANK. (B) SECTION A-A. 

FIGURE 20. - NUMERICAL MODELING OF CENTAUR LIQUID-OXYGEN PROPELLANT MOTION DURING SEPARATION MANEUVER. TIME, 1.0 SEC 
AFTER IMPULSE INITIATION; MAXIMUM VELOCITY * 50 cm/sec. 
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(A) CROSS-SECTIONAL VIEW OF LIQUID-OXYGEN TANK. (B) SECTION A-A. 

FIGURE 21. - NUMERICAL MODELING OF CENTAUR LIQUID-OXYGEN PROPELLANT MOTION DURING SEPARATION MANEUVER. 
TIME, 10 sec AFTER IMPULSE INITIATION: MAXIUMUM VELOCITY =12 cm/sec. 



(A) CROSS-SECTIONAL VIEW OF LIQUID-OXYGEN TANK. (B) SECTION A-A. 

FIGURE 22. - NUMERICAL MODELING OF CENTAUR LIQUID-OXYGEN PROPELLANT MOTION DURING SEPARATION MANEUVER. TIME, 
20 sec AFTER IMPULSE INITIATION; MAXIMUM VELOCITY ~ 12 cm/sec. 
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